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Abstract 



Dark matter may be composed of superWIMPs, superweakly-interacting massive particles pro- 
duced in the late decays of other particles. We focus on the case of gravitinos produced in the late 
decays of sleptons or sneutrinos and assume they are produced in sufficient numbers to constitute 
all of non-baryonic dark matter. At leading order, these late decays are two-body and the accompa- 
nying energy is electromagnetic. For natural weak-scale parameters, these decays have been shown 
to satisfy bounds from Big Bang nucleosynthesis and the cosmic microwave background. However, 
sleptons and sneutrinos may also decay to three-body final states, producing hadronic energy, 
which is subject to even more stringent nucleosynthesis bounds. We determine the three-body 
branching fractions and the resulting hadronic energy release. We find that superWIMP gravitino 
dark matter is viable and determine the gravitino and slepton/sneutrino masses preferred by this 
solution to the dark matter problem. In passing, we note that hadronic constraints disfavor the 
possibility of superWIMPs produced by neutralino decays unless the neutralino is photino-like. 

PACS numbers: 95.35.+d, 98.80.Cq, 26.35. +c, 98.80.Es 
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I. INTRODUCTION 



SuperWIMPs, superweakly interacting massive particles produced in the late decays of 
weakly-interacting massive particles (WIMPs), are promising non-baryonic dark matter can- 
didates l,y|- Well-motivated superWIMP candidates are the gravitino in supersymmetric 
models 1, M S B S and the first excited graviton in universal extra dimension mod- 
els iig. 

The supersymmetric possibility is realized naturally in supergravity with a gravitino 
lightest supersymmetric particle (LSP) and a slepton or sneutrino next-to-lightest super- 
symmetric particle (NLSP). (Throughout this paper, "slepton" denotes a charged slepton.) 
Both the gravitino and the NLSP have weak-scale masses. The NLSP freezes out as usual 
with a relic density near the observed value. However, after time t ~ 10^ — 10^ s, it decays 
through 

I IG 

i) ^ uG . (1) 

The gravitino then inherits much of the relic density of the slepton or sneutrino [s*], and its 
relic density is naturally of the right magnitude without the introduction of new scales. This 
is in contrast to other gravitino dark matter scenarios, where the gravitino is a thermal relic 
and the desired density is obtained by an appropriately chosen gravitino mass uIq ~ keV jof , 
or the gravitino is produced during reheating , where the desired dark matter density 

is obtained only for a fine-tuned reheat temperature Tr ~ 10^'' GeV. 

The decays of Eq. occur well after Big Bang nucleosynthesis (BBN). An immediate 
concern, therefore, is that they might destroy the successful light element abundance predic- 
tions of BBN. In fact, BBN is not the only worry — the Planckian spectrum of the cosmic 
microwave background (CMB), the diffuse photon background, and bounds on late time 
entropy production from the coincidence between BBN and CMB baryometry all impose 
significant constraints. As demonstrated in Refs. 0, 01 , however, these constraints on the 
electromagnetic (EM) energy released in the decays of Eq. exclude some of the weak 
scale parameter space, but leave much of it intact. Of course, at the border between the 
excluded and allowed regions, slight deviations from standard cosmological predictions are 
expected, providing possible signals in future observations. In fact, the existing anomaly in 
the ^Li abundance prediction of standard BBN may already be interpreted as a signal of 
superWIMP dark matter [3]. Such signals are particularly welcome, since superWIMP dark 
matter is so weakly interacting that it is impossible to detect through conventional direct 
and indirect dark matter searches. 

The previous work, however, neglected hadronic energy produced in WIMP decays. This 
was natural, as the WIMP decays of Eq. ((T)) contribute only to electromagnetic energy, as 
we will discuss below. However, the three-body decays 

/ ^ IZG , vWG 

V vZG , IWG , (2) 
produce hadronic energy when the Z or W decays hadronically. Hadronic energy re- 
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lease is very severely constrained by the observed primordial light element abundances 
E^ . [13L [14 Eij l. and so even subdominant hadronic decays could, in principle, provide 
stringent constraints. These three-body decays may be kinematically suppressed when 
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mi- — itlq < mw,rnz, but even in this case, four-body decays, such as / — > l^G —>■ IqqG, 
contribute to hadronic cascades and may be important. In this paper, we determine the 
hadronic branching fractions and compare them to BBN constraints on hadronic energy re- 
lease. Although we focus on the supersymmetric case, our results may be extended straight- 
forwardly with minor numerical modifications to the case of graviton super WIMPs in extra 
dimensions. 

In evaluating the constraints, there are two possible approaches. As the superWIMP relic 
abundance is automatically in the right range, a natural assumption is that superWIMP 
gravitinos make up all of the non-baryonic dark matter, with 

f^swiMP ^ 0.23 . (3) 

This is the approach taken here. Note that f^swiMP need not be identical to Qq — not all 
relic gravitinos need be produced by NLSP decays. Our assumption therefore requires that 
the other sources of gravitinos be insignificant. For the supersymmetric case, this typically 
imphes an upper bound on reheat ternperatures of ^ 10^° GeV Higher reheat 

temperatures may also be allowed [3, and are desirable, for example, to accommodate 
leptogenesis. For the universal extra dimension case, the requirement of insignificant Kaluza- 
Klein graviton production during reheating implies Tr < 1 to 10^ TeV, depending, in part, 
on the number of extra dimensions 

On the other hand, the Universe has proven to be remarkably baroque, and there is no 
guarantee that dark matter is composed of only one component. One might therefore relax 
the constraint of Eq. © on the superWIMP energy density and assume, for example, that the 
NLSP freezes out with its thermal relic density ^^nlsp- The superWIMP gravitino density is 
then l^swiMP = ("^g'/'"^nlsp)^nlsp- ^^i^ approach, the gravitino density may be low and 
even insignficant cosmologically. From a particle physics viewpoint, however, the viability 
of the gravitino LSP scenario is still worth investigation, as it has strong implications for the 
superpartner spectrum and collider signatures, independent of its cosmological importance. 

These approaches differ significantly, not only in the viewpoint taken, but also in their im- 
plications. Suppose, for example, that the gravitino and NLSP masses are both parametrized 
by a mass scale m^vsY- The NLSP number density then scales as 1/msusY if one as- 
sumes Eq. 0, but scales as msusY if one assumes a thermal relic NLSP density, since 
^NLSP ^ {^'^)~^ ^ ^susY) where {av) is the thermally-averaged NLSP annihilation cross 
section. Low masses are excluded in the former case, while high masses are disfavored in 
the latter. We consider both approaches to be worthwhile, but defer discussion of the ther- 
mal relic density approach, along with its very different consequences for the superpartner 
spectrum and implications for collider physics, to a separate study jlS^ . 

Assuming a fixed f^swiMP here, we find that gravitino superWIMPs provide a viable 
solution to the dark matter problem. We determine the allowed masses for the NLSP and 
gravitino. At the same time, we find that the hadronic energy constraint is significant 
and does in fact provide the leading constraint in natural regions of parameter space. We 
conclude that the analysis below must be done in any scenario with similar late decays. Given 
any standard cosmology, a significant thermal relic abundance of NLSPs will be generated, 
these NLSPs will eventually decay to superWIMPs, and the hadronic constraints we discuss 
below must be analyzed before the scenario may be considered viable. 

In passing, we note that another possibility is that the NLSP is the lightest neutralino. 
For a general neutralino, the decays x ^G, hG produce large hadronic energy release. 
These two-body modes are absent for photino-like neutralinos, where the only two-body 
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decay is 7 ^ 76*. However, even injthis case, one cannot avoid hadronic activity from decay 
through a virtual photon, 7 7*6* Q^G-, with a branching fraction of (9(10~^). Given 
the stringency of constraints on hadronic energy release, the neutralino NLSP case is highly 
constrained. Assuming that ^swimp makes up a sig nificant fraction of the dark matter, 
natural regions of parameter space are excluded [18|. In the current paper, we therefore 
focus on the more promising possibilities of slepton and sneutrino NLSPs. 

The paper is organized as follows. In Sec. |nj we summarize what is known about con- 
straints on EM and hadronic energy release in the early Universe. In Sec. IIIIl we give a 
detailed description of the gravitino superWIMP scenario with a slepton or sneutrino NLSP. 
In Sec. IIVI we determine the hadronic energy release from slepton and sneutrino decays and 
find the allowed mass regions. We summarize in Sec. |3 The relevant gravitino interac- 
tions, decay amplitudes, and three-body hadronic decay width formulae are collected in the 
Appendix. 



II. LATE TIME ENERGY INJECTION AND BBN 

The overall success of standard BBN places severe constraints on energy produce d by 
particles decaying after BBN. Given the baryon density specified by CMB measurements [l9|, 

This 
current 
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the prediction of standard BBN [20|] for deuterium agrees well with observations 
concordance is less perfect, but still reasonable, for ^He jl^,!!^!- On the other hand 
observations of ''Li 0, HE 11] are consistently lower than predictions; this is the leading 
anomaly in light element abundances at present. The abundances of ^He and ^Li also provide 
constraints, although these are subject to more theoretical uncertainties as discussed below. 

In this section, we review the constraints from these element abundances on energy release 
in the early Universe. Although a unified analysis is necessary for detailed conclusions, we 
review the constraints on EM and hadronic energy release in turn. A schematic summary 
of the discussion in this section is given in Fig. ^ 

The constraints on EM energy release are well understood, and we begin with a brief sum- 
mary of the results. (For details, see Ref. [l^.) High energy photons and leptons produced 
by late-decaying particles are quickly thermalized through interactions with background 
(BG) particles. High energy photons, for example, initiate EM cascades through processes 
like 77bg e~^e~ and •ycBG 7e, and the initial energy is quickly converted to very soft 
photons. The same is true for high energy electrons and muons. High energy taus decay 
to mesons, but these decay and produce EM cascades before interacting hadronically. (See 
Ref. 01 and the analysis below.) Injected leptons of all three families therefore contribute 
dominantly to EM energy and negligibly to hadronic energy. 

The energy distribution of the resulting photons is highly suppressed for > -Emax, 
where 



E„ 



22 T 



12 MeV 



keV 



(4) 



is related to the energy threshold for 77BG ~^ e"'"e 2^. Here T is the BG photon temper- 
ature; in the radiation-dominated era, it is related to time by 



t ~ 1.0 X 10*^ s 



keV' 



(5) 



4 



Time 


(sec) 1 10 102 103 104 105 106 10^ 


Hadronic Energy 


Thermalizes through 
EM interactions 


Interacts with background hadrons 


p + e n + v 
No constraints 


n+ p^D^ 4|-|g 
^He overproduction 


n+ p^D n +4He D 

D overproduction 


ElVI Energy 


Thermalizes through EM interactions 


Constraints weak 


y + D ^ n + p 
D overdestruction 


Y + ^He^D + D 
D overproduction 


Overall 


No constraints 


"He overproduction 


D overproduction 


Possible EM/had 
cancellation 


D overproduction 



FIG. 1: Summary of the leading constraints from D and ^He on late time EM and hadronic energy 
injection into the early Universe. Constraints from ^He and ^Li may also be important. (See text.) 



Because of this rapid thermalization, constraints on EM energy depend on the total amount 
of energy released and are essentially independent of the shape of the initial energy spectrum. 

For t ^ 10^ s, the BBN constraints are rather weak since -Emax is too small to destroy any 
light elements. For 10'* s ^ t ^ 10^ s, the Universe is hot and -Emax is low, and so only light 
elements with very low binding energies can be destroyed. The main constraint is from the 
overdestruction of D, with a binding energy of only 2.2 MeV. For t ^ 10''' s, the Universe is 
cooler and -Emax becomes high enough that ^He, with binding energy 19.8 MeV, is destroyed. 
Even a small amount of ^He photo-dissociation becomes a significant production mechanism 
for ^He and D, since the observed abundances of ^He and D are much lower than that of 
^He. The main constraint is from the overproduction of D. 

Between 10^ s and 10^ s there is a region where the overdestruction and overproduction 
of D cancel. For sufficiently high energy release in this region, the D abundance is that of 
standard BBN, but ''Li, with binding en ergy 2.5 MeV, is destroyed through ^Li-l-7 n+^Li. 



3rgy . 

The observed ^Li abundance [2^ |25|, 2y| is in fact low compared to the predictions of 
standard BBN 0. We see, then, that the ^Li abundance can be brought into agreement 



with observations without upsetting the D concordance for decay times in this window. For 



r] = UB/n^ = 6-0 X 10 the best fit values are j27 

^EM ~ 10"^ GeV , t ^ 3 X lOS , (6) 

where 

^EM = cem-Bem^nlsp (7) 

is the emitted EM energy density. Here eEM is the initial EM energy released in NLSP decay, 
and Bem is the branching fraction of NLSP decay into EM components. Vnlsp = ''^nlsp/''^7 
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is the NLSP number density just before NLSP decay, normalized to the BG photon number 
density = 2({3)T^/7t'^. It can be expressed in terms of the superWIMP abundance: 



Fnlsp - 3.0 X 10 



-12 



TeV 



^SWIMP 



SWIMP 



0.23 



(8) 



The parameters of Eq. © are naturally realized in the superWIMP scenario when a weak- 
scale NLSP decays to a weak-scale superWIMP with the observed dark matter density 0. 

The numerical constraints on EM energy release as a function of the particle decay lifetime 
have been determined in Ref. ji^], using the following criteria: 



D low 
D high 
^He low 
^Li low 



D/H < 1.3 X 10"^ 
D/H > 5.3 X 10~^ 
Yp < 0.227 
^Li/H < 0.9 X 10 



-10 



(9) 
(10) 

(11) 
(12) 



We adopt these constraints in our consideration of EM energy release bounds. 

Hadronic energy release, and its effects on BBN, has also been considered |ll|. El El H 
15l |. Although we are most interested in energy release at times t ^ 10"^ s, we summarize 
here what is known for t ^ 1 s for completeness. 

As we will show below, for the decay times of interest, mesons decay quickly and so 
contribute to EM energy. Only protons and neutrons are sufficiently long-lived that they 
might be able to interact hadronically. For times t ^ 150 s, however, even high energy 
protons and neutrons (with energies ^100 GeV) interact through EM scattering with BG 
photons before interacting with BG baryons. Injected high energy baryons then quickly 
become slow baryons, with the initial energy dissipated to EM energy. For such early times, 
this EM energy is harmless, as discussed above. Constraints are therefore functions of the 
number density of injected baryons -Bhad^NLSP, and independent of the initial baryon energy 

Chad. 

Slow baryons modify the number of neutrons and protons and therefore BBN predictions. 
In the very earliest era with t ^ 1 s, however, the weak interaction p + e^n + u is efficient, 
and so any deviation in the n/p ratio is washed out, and there are no strong bounds. For 
t > 1 s, the weak interactions decouple. Slow baryons produced in late decays cause the n/p 
ratio to deviate from the standard BBN value. In particular, extra neutrons create D. For 
1 s ^ t ^ 100 s, the resulting D is immediately burned to form ^He, and so the dominant 
constraint is from overproduction of '^He. This constraint is weak since the ^He abundance 
of standard BBN is already large. For t > 100 s, D can no longer burn into "^He 2| and 
the dominant constraint comes from D overproduction. This bound is significantly stronger, 
since the D/H ratio is ~ 10~^. 

At times t ^ 150 s, injected baryons interact with BG baryons before thermalizing. 
These interactions create secondary hadrons that can modify BBN predictions. In this era, 
constraints depend in principle on both ehad and -Bhad^NLSP separately. In principle, the 
constraints do not simply depend on the product, because the total number of secondary 
hadrons created by an injected hadron is not always proportional to the injected hadronic 
energy. The dominant constraint is from secondary baryons destroying ^He to overproduce 
D. 

Quantitative analyses on the BBN constraints of the hadronic energy injection have been 
performed in Refs. [ll, l^, IE 3 l3- discussed above, formally, the constraints depend 
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on the lifetime r of the decaying particle and on ehad and -Bhad^NLSP separately. However, the 
most recent analysis finds that for decay times r ^ 10^ s, when the hadronic constraint 
is most important, the constraints depend on the product ehad-Bhad^NLSP to within a factor 
of 2 for 100 GeV < Chad < 1 TeV We will therefore assume this behavior and apply the 
hadronic constraint to times r ^ 10^ s, where it is stronger than the EM constraint in the 
superWIMP scenario. 

The authors of Ref. assume the following constraints on the light element abundances 
(95% CL): 

D/H = (2.8±0.8) X 10-^ (13) 
= 0.238 ±0.004 ±0.010 (14) 
logio(^Li/H) = -9.66 ± 0.112 ± 0.6 . (15) 

In our analysis of hadronic energy injection, we take the constraints of Ref. derived 
using these bounds and 5had = 1- 

Note that the constraints of Eqs. ()13 |) - p5|) adopted for the hadronic energy analysis are 
less conservative than the constraints of Eqs. (j^-()12|) used in the EM analysis of Ref. j27l |. 
In particular, a significantly more aggressive D bound is assumed in the hadronic analysis. 
There is now impressive concordance between the baryon number determinations from D 
and CMB measurements, which further supports the narrow range of D/H given in Eq. (|13|). 
At the same time, existing anomalies in standard BBN may be indications that systematic 
errors are underestimated or that there is new physics, which would likely distort the D 
abundance, since D is very weakly bound. 

There are several further comments to make. First, we have been treating the EM 
and hadronic effects separately. Of course, these effects can add constructively or destruc- 
tively jl^]. We have explained above that the dominant hadronic constraint for t ^ 10^ s 
is from D overproduction, whereas for the EM constraint it is from D overproduction for 
t ^ 10® s and from D overdestruction for 10^ s ^ t ^ 10® s. For phenomena contributing to 
both EM and hadronic energy, then, the effects add constructively for t > 10® s, but may 
cancel for 10^ s ^ t ^ 10® s, and there may be fingers of allowed region in the latter time 
interval. 

Fina,lly, although we consider only the constraints above on D/H, Yp and ^Li/H, in 
Ref. constraints from ^He/D and ®Li/H are also analyzed for late time EM and hadronic 
injections. These two abundances, however, are subject to more uncertainties. ^He/D is 
considered because the destruction of '^He could produce ^He more easily than D, while the 
bindin g en ergy of ^He is higher than that of D. However, the interpretation of ^He measure- 
ments 3^1 may be subject to ambiguities arising from stellar production and destruction (3ll |. 



The ratio of ®Li/^Li has been observed in low metallicity halo stars and might also constrain 
late decays. The implications for EM cascades have also been explored in Refs. j3^ . l33l |. 
Again, however, there are arguably significant systematic uncertainties arising from the 
higher depletion of ®Li relative to ^Li as a result of its smaller binding energy j3J]. A better 
understanding of the connection between the primordial and observed abundances of ®Li 
and ^Li could provide a very powerful constraint on late decaying particles j3^. 

III. SLEPTON AND SNEUTRINO NLSP DECAYS 

In the superWIMP dark matter scenario, super WIMPs have no effect on the early thermal 
history of the Universe. They are produced only through metastable NLSP decays, NLSP— ^ 
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G + SM particle, giving a gravitino relic density of 



SWIMP 



"^NLSP 



th 

NLSP 



(16) 



^NLSP > 0-23 is allowed, since mg/mNLSP < 1- The SM particles produced in the NLSP 
decay induce EM or hadronic cascades, which will affect the primordial abundance of light 
elements, as discussed in Sec. |TT| 

The dominant decay mode for slepton NLSPs is / IG. Most of the energy released 
from the decay is transferred to EM cascades. The total energy released in slepton decay is 



m 



total 



NLSP 



rri; 



2mNLSP 

The total injected EM energy density is, for Am = mj 



(17) 



^EM = eEuBEuYj ~ 0(10-9 GeV) 



EM 



total 



Am 



300 GeV 



TeV 



m^ J 



SWIMP 



0.23 



(18) 



where eEM-^EM depends on the slepton's flavor, and Eq. (jH} has been used. For selectrons 
and smuons, energetic electrons and muons quickly thermalize with the background photons. 
All the visible energy is transferred to the resulting EM cascade: cem = -^totab -Bem = 1. In 
the case of staus, the decays of taus could also produce mesons and neutrinos, in addition 
to electrons and muons. The EM activity depends on the lifetime and energy distribution 
of decay products. As estimated in Ref. the EM energy released in stau decay is in the 
range 0.3Etotai to -Etotai- 

For the sneutrino NLSP case, the sneutrino mainly decays into neutrino and gravitino. 
The emitted high energy neutrino interacts with BG neutrinos and loses energy. If the 
neutrinos annihilate through uubg ~^ e~^e~, an EM cascade is initiated. However, the 
annihilation takes place through weak interactions, and so its cross section may be small. The 
effect of the expansion of the Universe must also be taken into account |3^,37j. Sneutrino 
decays induce far less EM energy release than sleptons, but for decay times r ^ 10^ s, there 
may still be some impact on BBN. For such decays times, however, much of the parameter 
space is excluded by considerations of hadronic BBN constraints, which will be discussed 
below. In what is left of this region, the gravitino mass is rriQ < (9(10) GeV. Such light 
gravitinos are less motivated in supergravity and so we have not considered this scenario in 
great detail here. 

The decay width of the slepton into LSP gravitino is 



r(/ 



L,R 



Il,rG) 



nir 



96nMi m% 



1 

2 



(19) 



where M4 = (IGvrGjv) For mQ/mj 1, the slepton decay lifetime is 

■100 GeVl^ ^ 



rilL,R ^ Il,rG) - 3.6 X 10« s 



Am 



m 



G 



TeV 



(20) 



This expression is valid only when the gravitino and slepton are nearly degenerate, but we 
present it here as a useful guide for the reader. In the analysis below, we always use Eq. f|T^ 



8 




1 5 10 50 100 5001000 1 5 10 50 100 5001000 



mg (GeV) mg (GeV) 

FIG. 2: (a) Lifetime (in seconds) for Il,r, v NLSPs and (b) EM energy release (in GeV) for the 
tr NLSP as functions of gravitino mass and mass difference 5m = mi — — mz- For (b), 
we take eEM = O.SE^totai, -Bem = 1, and Oswimp = 0.23. 



for the evaluation of the slepton decay hfetime. The slepton hfetime is given in Fig. in 

the {rriQ^ 5m) plane, where 

5m = Am — mz = mj — m^ — mz ■ (21) 

Note that in the definition of 5m, an additional mz is subtracted from Am; three-body 
decays, to be discussed below, are therefore always kinematically possible for 5m > 0. The 
typical decay time is 10^ s < r < 10^ s. For such decay times, there are strong constraints 
from the primordial abundance of the light element, as discussed earlier. 

Contours of EM energy release ^em = cem-Bem^" are given in Fig. ^ for NLSPs, 
assuming eEM = 0.5-Etotah -Bem = 1, and ^swimp = 0.23. We see that the EM energy release 
varies from 10^^ — 10~^ GeV and is largest for large 5m and small m^, when the decay 
lifetime is the shortest. The BBN constraints, however, are much weaker for early decays. 
As discussed in Ref. considering EM energy constraints only, one finds that for a stau 
NLSP, ~ 300 - 500 GeV and 5m ~ 200 - 1000 GeV are allowed. Moreover, late time 
EM energy release from stau decays to gravitinos may resolve the current discrepancy in the 
^Li for [m^, 5m) ^ (450 GeV, 200 GeV). For selectron or smuon NLSPs, the EM constraints 
are slightly stronger since all the released decay energy is transferred to EM cascades. 

The EM BBN constraints may be weakened by a number of possibilities. First, for 
decays in the range 10^ s < r < 10^ s, the EM constraints may be weakened by cancelation 
between the effects of EM and hadronic energy release, as discussed above. Second, if not 
all of non-baryonic dark matter is superWIMPs, I^swimp < 0.23, the contour lines for ^em 
will shift toward the left, since ^em oc I^swimp, and the allowed parameter space will be 
enlarged. Finally, if the NLSP is not a slepton but a sneutrino, the EM BBN constraints 
are, of course, much weaker, as the two-body decays contribute much less visible energy 
than in the slepton NLSP case. 
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FIG. 3: Feynman diagrams for slepton decays / I'VG leading to hadronic energy, x ™- diagram 
(c) is a neutralino or chargino. When kinematically accessible, V is a Z or W boson. For small 
Am < mzjmw, the leading hadronic decay is I — > lj*G, followed by 7* qq. 



IV. CONSTRAINTS FROM HADRONIC ENERGY RELEASE 

As discussed in SeclHl the abundances of Ught elements constrain the amount of hadronic 
energy injection into the early Universe. Although, for leptonic NLSPs, the hadronic energy 
release is much smaller than the EM energy release, the constraints on the hadronic energy 
release are much stronger. It is therefore important to evaluate in detail the hadronic activity 
induced by slepton and sneutrino NLSP decays. 

Let_ us first study the slepton NLSP possibility. The main decay mode of sleptons is 
I —>■ IG. For selectron and smuon NLSPs, such decays induce only EM cascades. For stau 
NLSPs, however, the resulting tau lepton may decay into mesons, and these might interact 
with background hadrons and induce dangerous hadronic cascades. The interaction time of 
the produced mesons is 



strong-int 



[((Xhad^^)^B] ^ ~ 18 S 



100 mb 



(0"hadf) 

1.8 X 10"^ s 



6 X 10 

7] 

100 mb 

(cThad^^) 



-10 



keV 



T 

6 X 10-1° 

7] 



1 s 



3/2 



(22) 



where (o"had'f ) is the thermally averaged strong interaction cross section. For the decay times 
of greatest interest here, 10^ s ^ r ^ 10^ s, this interaction time is long compared to the 



typical meson (vr,^^) lifetime, (-EMeson/'^^Meson) x 10~^ s. The produced mesons thus decay 
before they interact with background hadrons. 

As a result, the hadronic activity from the dominant two-body decay of slepton is neg- 
ligibly srnall. The main contributions to hadronic activity therefore come from three-body 
decays, / IZG, viWG with Z or W decaying hadronically, or from four-body decays 
/ l'y*G with 7* qq. The corresponding Feynman diagrams are shown in Fig. El For- 
mulae for the decay amplitudes and decay width calculations are listed in the Appendix. 

The hadronic branching fraction is defined as follow: 



B 



had 



T(J ^IZG)B^ + T{1 



l'WG)B^ + V{1 



I'qqG) 



V{1 IG) 



(23) 



where B^,B^ ~ 0.7 are the Z and W hadronic branching fractions. For right-handed 
sleptons, r{l —>■ uWG) = 0, and there are additional suppressions because the Z couplings 
are smaller than the W couplings. Thus, the hadronic branching fraction of right-handed 
sleptons is smaller than that of left-handed ones, especially when mw < Am < niz- 
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FIG. 4: Hadronic branching fractions for (left) NLSPs and (right) D NLSPs as functions of 
the NLSP decay lifetime. Along each curve, the NLSP mass is fixed at the value indicated and 
the gravitino mass varies. The neutralino/chargino parameters are chosen to be Mi = 2mNLSP) 
M2 = ^ = 4mNLSP) and tan/3 = 10. For the tr case, the dashed line gives the hadronic branching 
fraction resulting from tG, where the tau decays to a meson, and the meson interacts 

hadronically before decaying. We take the meson rest lifetime to be 3 x 10~^ s. 



For small Am, the slepton lifetime increases dramatically, since it grows as (Am)"^ 
by Eq. (pUj) . For Am < m^y, three-body decay modes are kinematically suppressed and 
four-body decay through a virtual photon dominates. The branching fraction of four-body 
decay is roughly expected to be (9((a;/47r)^) ~ O(10~^), which sets a lower limit for the 
hadronic branching fraction. Adjusting this estimation by taking into account details of the 
kinematics has negligible impact to our results. 

For the sneutrino NLSP case, the constraints from bounds on EM energy release are weak, 
especially at late times r ^ 10^ s, as discussed in Sec. IIIII However, the hadronic energy 
release is of the same magnitude for sneutrino NLSPs and for left-handed slepton NLSPs. 
The major hadronic activity comes from three-body decays, u uZG, IWG. The branching 
fraction for the four-body decay v i"~f*G —>■ vqqG is also similar to the analogous decay in 
the slepton case, except that diagrams (a) and (b) in Fig. |21 are absent since the neutrino is 
neutral. Last, we note that in the special case where the Bino and Wino masses are nearly 
equal, diagram (c) of Fig. El is highly suppressed, since the gauge boson V is almost purely 
Z. 

In Fig. |3]we show the hadronic branching fractions for tr and v NLSPs for a particular 
set of values of Mi, M2, /i and tan [3. Varying these parameters typically does not change our 
results significantly. The reason is that these parameters only affect the neutralino/chargino 
appearing in diagram (c) of Fig. 01 The contribution from diagram (c) is negligible if the 
neutralino/chargino is not degenerate to within about 10% with the sleptons. In addition. 



11 




FIG. 5: Hadronic energy releases ^had = Chad^had^NLSP as a function of decay lifetime tnlsp for 
(left) ffi NLSPs and (right) D NLSPs. Along each curve, ttIq is held fixed and Am = ttj-nlsp — 
varies. Particular values of Am are marked by the symbols indicated. The gravitino mass is 
itIq = 10, 10^, 10'^, 10^ and 10^ GeV for the curves ordered from left to right by their value at 
Am = 100 GeV. Constraints from BBN with and without ^Li/H, are also shown; regions 
above the curves are disfavored. We have assumed r^swiMP = 0.23 and ehad = ^("Inlsp — ""^g)- 



diagram (c) does not interfere with the other diagrams if all couplings are real, which elim- 
inates the possibility of cancelation. Hadronic branching fractions for left-handed sleptons 
are similar to those shown for sneutrinos. The branching fractions are smaller for the right- 
handed slepton case for the reasons mentioned above. In both cases, the hadronic branching 
fractions drop quickly for lifetime r ^ 10^ s, which corresponds to small 6m. For the 
case, mesons produced in r decay may, in principle, induce showers before they decay, as 
discussed above. This contribution to hadronic energy is also shown in Fig. |3J As can be 
seen, it is important only at early times r ^ 10^ — 10^ s and may be safely ignored for the 
decay times of most interest here. 

Note that the hadronic branching fraction of Eq. (|23p is defined at the quark level. 
Hadronization effects are important — for example, as we have shown above, mesons pro- 
duced in Z and W decays do not contribute to hadronic cascades since they decay before 
they interact. Such effects have been included in deriving the hadronic BBN constraints 
that we use [l^ . 

Fig. El shows the hadronic energy injection ^had = Chad-Bhad^NLSP as a function of decay 
lifetime tnlsp for f/j and u NLSPs. Along each curve, is held fixed and Am = ttinlsp ~ 
rriQ varies. The curves are truncated at Am = 100 GeV because for smaller Am, three- 
body hadronic decays are highly suppressed, leading to negligible hadronic energy injection. 
Hadronic BBN constraints from Ref. are also shown, with and without ^Li; the regions 
above the curves are disfavored. It is clear that for larger m^, larger values of Am becomes 
viable. Constraints are stronger in the i) case than in the case because the hadronic 
branching fraction is larger for left-handed particles than for right-handed ones. Notice that 
we have fixed f^swiMP = 0.23 in these plots. The lines could be simply rescaled for other 
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FIG. 6: Allowed and disfavored regions of the {mQ,5m) plane, where 5m = w-nlsp 
for (left) and (right) D NLSPs. The shaded regions are disfavored by CMB, EM BBN, and 
hadronic BBN constraints, as indicated. The lines correspond to EM constraints from ^He/D 
(solid) and hadronic constraints from ^Li/H (dashed). The circle indicates the best fit region 
where the ^Li discrepancy is resolved. See text for details. We have assumed riswiMP = 0.23 and 
Chad = 5('TiNLSP ~ ''tt-g)- Note that the axes have different scales in the left and right panels. 



values of ^swimp since Inlsp oc I^swimp- We have also adopted the hadronic BBN analysis of 
Ref. If updated analyses of the BBN constraints become available, it is straightforward 
to impose it in these plots to find out the remaining viable parameter regions for itlq and 
Am. 

Fig. ini shows the allowed and disfavored regions of the {rriQ, 6m) plane for the two NLSP 
cases. The CMB distortion constraint, the EM BBN constraint of Ref. jl^], and the 
hadronic BBN constraint of Ref. are included. For the f/j case, the CMB constraint 
excludes a region with 6m < 150 — 300 GeV. EM BBN constraints from D /H+Yp disfavor 
regions of 6m below 400 GeV for m^ < 500 GeV and 6m below about 200 GeV for larger 
mQ. Hadronic BBN constraints from D/H+Yp disfavor the part of the parameter space with 
large 6m where the slepton decays relatively early. Note that the hadronic constraint is 
extremely important — it disfavors a large and natural part of parameter space that would 
otherwise be allowed. 

Even after all of these constraints, however, the unshaded area in the region m^ ^ 
200 GeV and 200 GeV ^ 6m ^ 1500 GeV still remains viable. The mass must be 
above 500 GeV, but this is within reach of the LHC. The "best fit" point where the ^Li 
discrepancy is fixed by EM energy release from late slepton decays is indicated by a circle 
at {m^, 6m) ^ (450 GeV, 200 GeV). Note that it also remains viable, even given the strong 
hadronic BBN constraints. 

We also indicate the region that is disfavored by the less solid, but stronger, constraints 
from ^He/D and ®Li/H. The ^He/D constraint excludes parameter space below the solid 
line. The allowed region would then be pushed up to larger mQ ^ 3500 GeV and larger 
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6m ^ 800 GeV. Such a heavy stau NLSP is beyond the reach of the next generation hadron 
colhders. On the other hand, taking the ^Li/H constraint hterally would exclude parameter 
space above the dashed line. Taking both constraints literally and combining them would 
exclude most of the region of parameter space shown. Clearly, a more firm understanding 
of BBN will have important consequences for the stau NLSP super WIMP scenario. 

We have assumed ^swimp = 0.23; of course, the allowed parameter space would be 
enlarged if I^swimp < 0.23. We have also taken ehad = II'^nlsp — i^g) ™- analysis; the 
allowed region could shift if there are significant corrections to this relation. 

The sneutrino NLSP case is also shown in Fig. IHl Hadronic BBN constraints from the 
D/H and Yp abundances only disfavor 5m ^ 100 GeV. Even including the stronger but more 
speculative ^Li/H constraint, there is still a large region of {mQ,6m) that is allowed. The 
origin of this result is that at r ^ 10^ s, the sneutrino NLSP is completely free of hadronic 
BBN constraints. We have again assumed I^swimp = 0.23. To get this relic density from 
a sneutrino thermal relic density, one would typically require sneutrino masses above 500 
GeV. However, even assuming a thermal relic abundance, if the superWIMP component is 
only part of the present dark matter, smaller sneutrino masses are allowed. 

For right-handed selectron and smuon NLSPs, the results will be similar to the right- 
handed stau NLSP case. The EM constraints are slightly stronger, since all the released 
energy induces EM cascades. For left-handed slepton NLSPs, the region of (m^., Sm) that is 
disfavored by hadronic injection is similar to the region disfavored for the sneutrino NLSP. In 
addition, however, there are EM BBN constraints that are similar to those for right-handed 
slepton NLSPs. 



V. CONCLUSIONS 

In this paper, we explored the possibility of superWIMPs as candidates for dark matter. 
Examples of superWIMPs are the gravitino in supersymmetric models and the lightest 
Kaluza-Klein graviton in universal extra dimension models. SuperWIMPs obtain their relic 
density through late decays: WIMP superWIMP -|- SM particles. The decay usually 
occurs between 10^ s and 10^ s, which has important cosmological implications. Such late 
decays release EM and hadronic energy into the Universe, which may affect BBN predictions 
for the primordial abundance of the light elements. The constraints on EM injection have 
been studied in detail in Ref. [2]. In this paper, we have analyzed the hadronic BBN 
constraints. We have taken the lightest gravitino as a concrete example of a superWIMP 
and have focused on slepton and sneutrino NLSPs as the most promising WIMP candidates. 

We have determined the hadronic energy release by calculating three-body decay widths 
in detail. For the cases of and u NLSPs, the hadronic decay branching fraction is below 
the level of 10"'^ when mf^ < 1 TeV and mc < 300 GeV, respectively, or when the decay 
lifetime is r ^ 10^ s. We identified the allowed and disfavored regions of the {mQ, 6m) plane, 
imposing CMB constraints, EM BBN constraints, and hadronic BBN constraints. For the 
sneutrino NLSP case, 6m ^ 60 — 200 GeV is allowed for a large range of mQ. For the stau 
NLSP, the hadronic constraints are weaker: 6m ^ 200 — 1200 GeV. However, additional 
constraints from CMB and EM energy injection apply. Combining all the constraints for 
the right-handed stau NLSP, the allowed window is 300 GeV < m^ < 1 TeV, 200 GeV < 
6m < 1200 GeV, corresponding to mf^ > 500 GeV. We have assumed ^swimp = 0.23 in our 
analysis. The constraints would be relaxed if superWIMPs are only part of the present dark 
matter. In addition, there are still ambiguities in the BBN constraints if the decay time lies 
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in the region where EM and hadronic effects are comparable. In particular, as we discussed, 
their effects on D might cancel, and the allowed parameter space could be enlarged. Progress 
in firming up BBN constraints in this region, and for other elements, such as ^He and ^Li 
may have crucial implications for the superWIMP dark matter scenario. 

Although the superWIMP itself would escape all direct and indirect dark matter searches, 
the slepton/sneutrino NLSP will have rich implications for collider phenomenology. Such 
metastable sleptons will not decay inside the detector, resulting in signals of highly ionizing 
tracks for sleptons and missing energy signals for sneutrinos. Discussion of the collider phe- 
nomenology, combined with thermal relic density calculations for the NLSP in supergravity 
models will be presented in Ref. [Ts^ . 

Although we have focused on the SUSY scenario in this paper, superWIMPs could in 
general be any gravitationally interacting particle that obtains its relic density through the 
late decay of a WIMP. The discussion and results for late decay of a leptonic WIMP in 
other models will be qualitatively similar to the discussion and results presented here. This 
scenario also suggests that a particle that appears to play the role of dark matter at late 
times, even after BBN, could very well be different from the particle that constitutes dark 
matter now. The very late decay of "would be" dark matter particles may have important 
cosmological implications, for example, affecting small scale structure. This feature r nay 
provide qualitatively new possibilities for explaining puzzling cosmological observations [38|, 

We have assumed i?-parity conservation in our discussion. In the case of -R-parity vio- 
lation (RPV), the gravitino could still constitute dark matter as long as its decay lifetime 
is comparable to or longer than the age of Universe 41 1. Gravitino dark matter in an 



RPV scenario could be distinguished from the superWIMP scenario by both cosmological 
observations and collider experiments. For example, the decay of even a tiny amount of 
gravitino dark matter into SM particles in an RPV scenario could be seen in the diffuse 
photon flux. In addition, if the RPV is not extremely small, the collider signatures could 
be different from those in the superWIMP scenario. Further work is needed to study how 
to distinguish these two scenarios. 
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APPENDIX: WIDTHS FOR THREE-BODY SLEPTON DECAYS 
1. Interactions 

The gravitino-gaugino-gauge boson interaction is 

^ = -^?.[7^7l7'^AF.p, (A.l) 
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where = A/2M4 = 2.4 x 10^^ GeV is the reduced Planck mass, and ipfj_ = ip^{x) denotes 
the gravitino field at spacetime point x. The gravitino-slepton-lepton interaction is 

where subscript h = L,R denotes the lepton chirality. In addition, there are gravitino- 
slepton-lepton-gauge boson interactions: 

C = ^^AMYY^, . (A.3) 

Here gy is the gauge coupling coefficient, which is given explicitly below. 
Since gravitinos are Majorana particles, 

= C^l = . (A.4) 
The gravitino spin sum (with four momentum p^) is 

T^GUpMp) = - h-p + ma) L.-^) 

h \ "^G / 



2. Matrix Amplitudes 

Here we present_the matrix amplitudes for the three-body decays / — > I'VG, where I' is 
the SM partner of / or its weak isospin partner, and V = Z or W. The Feynman diagrams 
are shown in Fig. El 

We define the following ten operators: 





= I', pr G'^ p, ■ e* 


(A.6) 


02 


= V^pi-G^pi-e* 


(A.7) 


03 


= I'hPc- Pi ■ e* 


(A.8) 


0, 


= lYy^ ■ e* pv ■ G" 


(A.9) 


05 


= T^^.py ^-e* pi-G'' 


(A.IO) 


0, 


= il'h 1 ■ Pv 1 ■ (^* Pv ■ G" 


(A.ll) 


O-j 


= il'hl-PGl-^* Pv-G" 


(A.12) 


Os 


= i\e* -G" 


(A.13) 


O9 


= i\^-py e*-& 


(A. 14) 




= ^ 7 ■ 7 ■ Pv-e* ■ G'^ , 


(A.15) 



where e* is the polarization of gauge boson V . Notice that for an on-shell gravitino in the 
final state, (^3 = by using the gravitino equation of motion. 
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a 



In terms of these operators, the matrix elements for three-body slepton decay [diagrams 
)-(d) in Fig. E] are 

(A.16) 

+ (A. 17) 



(d) in Fig 


E] are 








Ml^ = 


ih 
egv 


1 


(20i 




V2M, 




M\ = 


ih 






1 


V2M, 


ml 


+ mf. 





'^23 



V — 



^ '^iMTm' -m^ - (^9) - 4 + O, - O,,) + mlO^] (A.18) 



Ml = ^^Os , (A.19) 

where Xi denotes a neutrahno or chargino, the couphngs gy, g^^^, and g^^v are given below 
for each decay process, and 

^12 = (Pg + Pi'f ' "^13 = iPc + Pvf , = {Pv + Pvf ■ (A.20) 

Notice that + '^is + "^23 = + + ^1- Given this relation, the four-momentum 
scalar products that appear below in the expressions for squared matrix elements can be 
expressed in terms of m^2) "^is ^i^d rn^^. 

For each specific decay, the couplings in Eqs. ()A.16|) - ()A.19|) should be replaced according 
to the following rules: 



IZG 

h ^ L , V ^ Z (A.21) 

g{^ - g'i = ■ n (l - sin' Qw\ (A.22) 

sm By/ cos b'vi/ V 2 / 

^ g^. = . / —m COS 9w + N*, sin 9w) (A.23) 

^ V 2 sm 61^/ cos 6w 

giv ^ ^7x.^ = iV* (-sin^H/) + Ar;^(cos^H^) (A.24) 
■ vWG and vl ^ /VTG 

L , 1/ ^ (A.25) 
\/2 sin 



- ^7^^ = - . ^ (A.26) 



> ^ g'^_ = —1 V*^ (A.27) 

^» V2sin^vi/ 



h ^ R, V ^ Z (A.29) 



sm Uw cos O'ly ^ ^ 

9x. ^ ^?x" = ^^^1 (^-31) 

-^^ ^» cos Ow 

9x.v ^ 4^ = A^^i(-sin^H^) + Ar,2(cos^H') (A.32) 
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V (A.33) 



z/L vZG 



sm b'vi/ cos ow \ 2,/ 

g'i^ ^ g^^^ = ^ . / — (iV;^cos^H/ + iVasinM (A.35) 

^» V 2 sm 6w cos fclvi/ 

^ = iV* (-sin^H/) +iV;^(cos^H/) (A.36) 

In these expressions, V and are matrices that diagonahze the chargino and neutrahno 
mass matrices, following the conventions of Ref. ji^ . 



3. Squared Matrix Elements 

The differential decay width is 

dT 



(27r)3 32mj 



(A.37) 



where M = M^y + M\r + M^y + Mir. 

The sum of the matrix elements can be written as 

M{lh I'VG) =Ml^ + M\, ^Wy + M'^ 



i=1...10 



(A.38) 



where the M.i can be read off from Eqs. ()A.16|) - ()A.19|) above. The squared matrix element 



IS 



Mijh I'VG) 



i=1...10 

i<j i<j 

+ E ^eiM,M*)Ol^j+ E "^rn {MiM;)0';] . {A.39) 

i,j=1...10 i,j=1...10 



In our calculations, we chose the convention that the diagonalizing matrices V and A^ are 
real. All couplings appearing in Eq. ()A.16|) — ()A.19|) are therefore real. Thus, all the Aii 
are real except for Aig, which has both real and imaginary components. The only non- 
zero contributions to the last term in Eq. ()A.39|) come from Qm (A^jA^g) {i < 8) and 
'^m (MsM*) (j > 8). Expressions for d^i, O^^j, Of^ {i < 8) and O'^j (j > 8) are: 



Sm'^my 



m 



G 



2,2 



X 



X 



m^my + mjmy - 2myp^ ■ pi - {p^ ■ py) + 2p^ ■ pypj ■ py - {pj ■ py) (A.40) 



Sm^rriy 



^g'Pg-Pi+Pg-Pv mjmy-{pi-py 



(Pg ■ Pif + '^Pg ■ PiPg ■ Pv - (Pg ■ Pvf + m% {^1 + ^v- '^Pi ■ Pv 



O: 



3,3 



(A.41) 
(A.42) 
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X 



m|my - rrivPa ■Pi + '^ (Pg ' Pvf +Pg-Pv (Sm^ - 2pi ■ pv) 



3m% 



X 



m%my - rriyPa • - 4 (p^ • pvf +Pg-Pv {-^rnl + ^Pi ' Pv) 

^ m%m\ - {pq ■ Pv f 



3m% 



X 



m%my - rUyPQ • - 4 (p^ • pvf +Pg-Pv (-3m^ + 4p;- • pv) 



m^my -{pg-PvY 

2 



[m%my + 2m| (p^ • pvf + ^Pq ■ Pv - Pg ' Pi [m%my + 4 (p^j • pvY 
-m%PG ■ Pv {my - 2pi ■ py) } 

m%-PG- Pi +Pg-Pv 2m%my + (p^ • Pv) 



^m^rriy 



Sm^rriy 



2m ^my + {PG-Pvy 



X 



m%my - rriyPG • Pf - 2 (p^ ■ Pvf - Pq ' Pv {ruy - 2pj ■ py) 
^2mgmy + {pG-pvf 



3m^my 



{m'^G^y - 2m| {pg ■ Pvf - 4 (p^ • pvf 

+Pg ■ P'l [-m%my + 4 (pg ■ py)^] + m|p^ ■ Pv {^my - 2pi ■ py^ } 
m%-PG-Pi+PG-Pv 



3m'^my 



X 
X 



(Pg ■ Pi) + Pg ■ P«Pg • Pv + "1(5 V"'^« ^ Pi ■ 
mfmy + ttiIpg ■ Pj - Pg ■ PvPi ■ Pv + (P[ • PvY 



0\\ = 



3m^ 



m|mf - {pg • Pif] [- (Pg • Py)^ + m| (ml - pf ■ py) 
+Pg ■ Pv {-mf + Pi ■ Pv) + Pg ■ Pi (-ml +PG-PV+ Pi ■ Pv) 

OTa = 

{{Pg ■ Plf [m%ml + 2 (m? - ml) p^ • py 

G 

+2 (Pg ■ ("^y - Pi ■ Pv) - '^Pg ■ Pi [(Pg ' Pvf (mf - pj ■ pv) 
+m| (m? - PI ■ Pv) (my - pi ■ pv) + m}^PG ■ PvPj ■ Pv 
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+m| (^-m\m\ + {pi ■ Pv)^)] } 
oi, = = Ol% = O 



2,10 











rrin 



G 



m\my - {pq ■ pvf niy+pQ-pv- Pi ■ Pv 



rriri 



G 



rriQ-PG-Pi+PG-Pv 



v-iPG-PvY 



8 



SruQUiy 



my I m 



"G 



2Pg ■ Pi) Pg-Pv+Pg- Pv 



+m%my (my - pi ■ pv) - (Pa ■ Pvf {ruy + pi ■ pv) 



Sm^rriy 



m^my 



+ m%myPQ ■ Pv + m|my (p^ ■ Pv f - 2 (p^ • PvY 



-myPa ■ pi (mym^y + (p^ • pvf) - [pq ■ Pvf {my - "^Pf ■ Pv) 

- {m^my - m| {p^ ■ pvf - Sp^ ■ Py + Pa ■ Pi [-m},my + 3 (p^ 



+m^^PQ ■ Pv (Smy - 2pi ■ py) } 

/-/ire 
^b,7 



/rjre 



/rjre 



3m% 



m 



G 



my - {Pc-Pvf 



3m| 



2>m% 



(3m| - 2pQ ■ p^ Pq-pv + 2 (p^ ■ pvf + m% (rriy - pi ■ pv 
^ ^ni\^m\ - m\ {Pc'Pvf + {Pg 'Pvf 
-Pg • Pi {m%my + {pq ■ Pvf) - m^pQ ■ pv (rriy - 2pi ■ py) 
(3m| - 2p^ ■ pi) Pg-Pv-^ {Pg " Pvf 
(Pg • Pvf {rriv - 3pf • Pv) + m'^^my (rriy - pi ■ pv) 
— ^^^-myp^ ■ Pv - (m| - 2pg ■ p;-) {p^ ■ Pvf - 2 (p^ ■ pvf 

+m%my (rriy - pi ■ pv) + m| (p^ ■ pvf (rriy - pf ■ pv) 
3^2^^2^ [^I^vPg ■ Pv - {ml - 2pc • Pi) {Pg ■ Pvf - 2 (pg • pv) 
+m%my (my - pi ■ pv) + m| (p^ • pvf (rriy - pf ■ pv) 
m%m\ - m\ (p^ • pvf + (Pc • Pvf 



3m, 



G 



-Pg • Pi {m%my + (p^ • pvf) - m^PG ■ Pv {rriy - 2pi ■ py) 
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Ol 



10 



/■r\im 

^6,8 



r64|44 ,42/ \5 

Im^my + m^myp^ ■ pv + m^my [p^ ■ pv) 



-Pg ■ Pi 



m^my + m^niy {p^ ■ pv) - 4 {p^ ■ pv) 



+^6 (Pg ■ Pvf (Sm^ - 2pi ■ Pv) } 



2m^my + {p^ ■ Pv) mv+Pc-Pv -PrPv 



3mQmv 



Sra^mv 



2m^my + {PG-Pvy 



(m| - 2p^ ■ pi) p^-pv-2 {pq ■ pvf + m| (niy - pj ■ pv) 



3m, 



G 



g-Pg-Pi + Pg- Pv 2m^my + {p^ ■ Pv) 



(A.71) 
(A.72) 

(A.73) 
(A. 74) 



3m^m^ 



g'"v 



(m| -Pg-Pi+Pg-Pv){ [mfmy - {pi ■ Pvf + [pq ■ Pv){Pi ■ Pv) 



-mlipG ■ Pi)^ + {Pg ■ Pi){.Pg ■ Pv)iPi -Pv-Pg- Pv)} 



Sm'^my 



TZZri^G - Pg ■ Pi + Pg ' Pv) {rriimy - {pi ■ pv) )m^ - myipc ■ Pi) 



-{Pg ■ Pv)^Pi ■ Pv + {Pg ■ Pi) (Pg ■ Pv){ml + Pi ' Pv) 



^4,8 



3ml 



r { (Pg ■ Pi - ^Dpi ■ Pvm% - Pg ■ PjiPd " Pv) 



+Pg ■ Pv [pg ■ Pi^d + (Pi ■ Pv - 2mf)ml + {p^ ■ p;-)^] } 
= Ofg = 0fio = • 



(A.75) 



(A.76) 
(A.77) 



(A.78) 
(A.79) 
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